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Abstract 

l-D systems with deterministic disorder, such as those with quasiperiodic or 
substitution sequence potential distributions, have been much studied. It was 
recently shown that a generalization of quantum Turing machines (QTMs), 
in which potentials are associated with elementary steps or transitions of the 
computation, generates potential distributions along computation paths of 
states in some basis B, that are computable and are thus periodic or have 
deterministic disorder. These generalized machines (GQTMs) can be used 
to investigate the effect of potentials in causing reflections and reducing the 
completion probability of computations. This paper expands on this work by 
determining the spectral and transmission properties of an example GQTM, 
which enumerates the integers in succession as binary strings. A potential is 
associated with just one type of step. For many computation paths the po- 
tential distributions are initial segments of a distribution that is quasiperiodic 
and corresponds to a substitution sequence. Thus the methods developed in 
the study of l-D systems can be used to calculate the energy band spectra 
and Landauer Resistance (LR). For energies below the barrier height, The LR 
fluctuates rapidly with momentum with minima close to or at band-gap edges. 
Also for several values of the parameters involved there is good transmission 
over some momentum regions. 
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I. INTRODUCTION 



The discovery of quasicrystals PL (See for recent reviews) has stimulated much study 
of systems with deterministic disorder jp. These systems are neither random nor periodic. 
As part of this work there has been much interest in the study, either with Kronig-Penney 
or tight binding Hamiltonians, of the spectra and transmission properties of 1-D systems 
whose potential distributions correspond to substitution sequences. (These are defined in 
§). For the Period Doubling, |§|||7|], Thue- Morse and Fibonacci § substitution 

sequences, the spectra are singular continuous and are Cantor sets. For the Rudin-Shapiro 
sequence, the spectral properties are unknown ||. Potential distributions corresponding to 
circle sequences |10| and the prime number distribution JTl[] have also been studied. 

Transmission properties of 1-D structures for which the potential distributions are initial 
segments of substitution sequences have also been studied |^2]-|16|. Of special interest here 
are studies of the Landauer Resistance JT7| , [T%| ] for various potential distributions and its 
relation to the band structure of infinite periodic systems whose unit cells are initial segments 



of substitution sequences [19-21 



One dimensional structures of these types also appear in a generalization of quantum 
Turing machines introduced elsewhere p2| . The generalization consists of modifying the 
Hamiltonian description of quantum Turing machines (QTMs) so that potentials are asso- 
ciated with some types of steps or transitions of the computation as it moves along paths of 
states in a suitable basis. The resulting distributions of potentials along these computation 
paths are computable and are thus either periodic or deterministically disordered. As a 
result generalized quantum Turing machines (GQTMs) are both quantum computers and 
1-D structures of the type of much recent interest. 

The presence of these potentials means that GQTMs also serve as a test bed to examine 
the effect of these potentials on the transmission and reflection of the computation state along 
any computation path and the completion probability for the computation. The importance 
of this for quantum computation has been repeatedly emphasized by Landauer |23| . He has 
noted that these potentials, which may result from environmental influences and errors in 
construction of actual physical models of the computation process |2"3| . |2"5| ] , cause reflections 
back along the computation path and may result in exponential decay or localization of the 
transmitted state. 

This emphasis is based on the fact that in much of the work done so far on quantum 
computation, the computation as been assumed to proceed smoothly along computation 
paths with no potentials present. This is the case whether quantum computers are used 



as computers |26|-j33|1 or as simulators of other quantum systems |3Jj3^]. Work on Shor's 
quantum algorithm jSBJ and on quantum gates and quantum error correction PT| , |55|] is also 
included. 

The main purpose of this work is the examination of the effect of potentials on the spectral 
and transmission properties of computation states along computation paths. The method 
will be based on other work in which potentials are associated with some types of steps 
or transitions of the computation as it evolves along paths of states in a suitable basis. The 
Feynman Hamiltonians for these GQTMs can, under certain conditions, be decomposed into 
a direct sum of tight binding Hamiltonians, H = J2i Hi, where the i sum is over computation 
paths of states in a basis. The potential distribution, Vi, in Hi depends on the path i. 
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Here the example GQTM, which was introduced in other work [Q, will be examined in 
detail. This example, the counting GQTM, which enumerates the integers in succession as 
binary strings, has a potential associated with read-1 steps only, (read-1 steps are defined 
at the end of Section |T|.) For many paths the potential distribution was seen to be qua- 
sicrystalline and to correspond to initial segments of a substitution sequence generated from 
an infinite alphabet. 

Both the spectral and transmission properties of the tight binding Hamiltonians asso- 
ciated with the counting GQTM will be determined. The energy band spectra will be 
calculated for associated periodic systems. These are 1-D crystals whose unit cells are (fi- 
nite sections of) computation paths that contain the potential distributions being studied. 
Transmission properties will be calculated by determining the Landauer Resistance of com- 
putation paths for the counting GQTM. This paper should be read in conjunction with |22[ 
as much of the earlier work will not be repeated. 

In the next section the physical model for QTMs (generalized or not) will be presented. 
This will be followed in Section [Ifl] by a brief description of the generalization of QTMs 
and the decomposition of Hamiltonians for GQTMs into a sum over paths of tight binding 
Hamiltonians. Computation paths will be defined. The counting GQTM will be described 
in Section [TV[ The discussion of the step operator and potential distribution will be brief 
as much of the material is already given elsewhere [22|. The transfer matrix will be defined 
for potential distributions containing the first 2 Ar_1 potentials with N an arbitrary positive 
integer. 

The emphasis of Sections [V| and [V| is on the spectral gaps and band widths for the 
associated periodic systems and Landauer Resistance of systems described by tight binding 
Hamiltonians with 2 N ~ l potentials for the counting GQTM. The width and separation of the 
potentials is related to an initial segment of a generalized substitution sequence described 
elsewhere p2| , ^9| . The work is exploratory with emphasis on calculations to show the de- 



pendence of these properties on various parameters. Also as was done by Roy and Khan 
|20|| , the relation between the Landauer Resistance of a system and the band spectra of the 
associated periodic system will be investigated. 

It should be noted that the only effect of this generalization is the introduction of poten- 
tials along computation paths. There is no change in the class of QTMs or in the concept 
of quantum computability. 



II. THE PHYSICAL MODEL 

The physical model for one-tape quantum Turing machines consists of a two-way infinite 
one dimensional lattice of systems each of which can assume states in a finite dimensional 
Hilbert space. If the space is two dimensional, the systems are referred to as qubits. This 
term will be used here even if the dimensionality is greater than two. It is often convenient 
but not necessary to consider the lattice as spin systems, e.g. spin 1/2 systems for binary 
qubits. 

A head which can assume any one of a finite number of orthogonal states |Z) with leL 
moves along the lattice interacting with qubits at or next to its location on the lattice. 
Elementary QTM actions include one or more of (1) head motion one lattice site to the 
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right or left, (2) change of the state of the qubit scanned or read by the head, (3) change of 
the head state. What happens depends on the states of the head and scanned qubit. 

Here the system states are all assumed to lie in a separable Hilbert space 7i. Based on 
the above description a particular basis, the computation basis, defined by the set of states 
{\l, j, S)} and which spans TC, is used. Here l,j refer to the internal state and lattice position 
of the head. The qubit lattice computation basis state \S) = ®m=-ao\S( Tn )) where S is a 
function from the integers to the qubit state labels (e.g. {0, 1}) such that S(m) ^ for 
at most a finite number of values of m. This condition, the tail state condition, is one 
of many that can be imposed to keep the basis denumerable. Note that all qubit states, 
such as a|0) + 6|1) for arbitrary a, b and sums of products of these states are included as 
appropriate linear superpositions of the computation basis states. 



III. TIGHT BINDING HAMILTON! ANS FROM QUANTUM TURING 

MACHINES 

One begins by associating a bounded linear operator T with each GQTM. T is referred 
to as a step operator for the GQTM because iteration of T (or its adjoint) corresponds to the 
successive steps in the forward (or backward) time direction of the GQTM. An infinitesimal 
time interval is associated with the steps of T so that it can be used directly to construct 
a time independent Hamiltonian. The specific construction used here is that proposed by 
Feynman |[40|| : 

H = K(2 - T - T f ) (1) 

where K is a constant. This definition has the advantage that if T = Y where Y is the 
bilateral shift along the lattice, then H is the kinetic energy of free head motion on the 
lattice. As such it is equivalent to the symmetrized discrete version of the second derivative, 
(-h 2 /2m)d 2 /dx 2 . 

For each GQTM the associated step operator T is defined as a finite sum over elementary 
step operators. That is T = J2i, s Ti s where the sum is over all head state labels leL and all 
qubit state labels seS where S is the set of qubit state labels. For binary qubits S = {0, 1}. 
Ti 7 g corresponds to the action taken by the GQTM associated with T when the head in state 
|/) sees or reads a qubit in state \s). 

The generalization consists in noting that for each /, s Ti >8 = 7z jS H / z, s where 7^ is a 
positive real constant < 1. For QTMs without the generalization 7^ = 1 for all I, s. Wi jS is 
defined as a sum over all lattice positions of products of projection operators, bit and head 
state change, and head position change operators. The specific definition of W\ s appears 



elsewhere 22 



From now on either the W\ a or the TJ jS will be referred to as the elementary step operators 
for the GQTM. The reference wiil be made explicit in contexts where the difference between 
Ti tS and W\ >a is important. 

The main condition imposed on T and its adjoint is that they be distinct path generating 
in some basis B which spans TC. This means that iterations of T or on any basis state 
in B, generate a path of states that, up to normalization, are also states in B. Furthermore 
the paths are distinct in that no two paths join, intersect, or branch. 
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Although B can be any basis which spans H, (see p7| for details), B is required here to 



be the computation basis, defined in Section [E|. A path in B is defined as a set of states in 
B ordered by iteration of T or If the state \l m , j m , S m ) = \m, i) is the mth state in some 
path i, then the m + lst and m — 1th states in the path are given respectively by \m + 1, i) = 
Wi m , Srn (j m )\l m Jm,S m ) = \l m+1 , j m+1 , S m+1 ) and \m-l,i) = W^^g^^^lm, j m , S m ) = 
\lm-h jm-h Sm-i)- These equations must be modified if either T or annihilate \m,i). 

Mathematically the requirement that T be distinct path generating in the computation 
basis, B, can be expressed by the condition that T is a direct sum of weighted shifts on 
B. That is, T = UD = Yl,iP%Pi where U = J2i Ui is a direct sum of shifts (bilateral, 
unilateral, adjoint of a unilateral shift, finite, or cyclic) and D is pure discrete with positive 
real eigenvalues and eigenstates in B. These types of shifts correspond respectively to the 
associated paths being two-way infinite, one-way infinite from either end (m > or m < 0) 
finite with distinct ends, or cyclic. One also has U = J2i, s Wi tS . 

For each i Tj = UD is a weighted shift on a subset Bi of B which spans a subspace 



PiH of H. pi|| . The projection operators Pj are all pairwise orthogonal. More explicitly if 

{\m,i)} is the basis Bi and Ui\m,i) ^ 0, Ti\m,i) = Di^ m \m + D im is the eigenvalue 
of D associated to the eigenstate \m,i). For the adjoint T^\m,i) = Dj jTO _i|m — if 
Uj\m,i)^0. 

These results, expressed in the computation basis labels, give for any state which is not 
a terminal path state, 

T\l m ,j m , S m ) = Dl m ,S m (jm) IWli 3m+l, Sm+l) 
T \lrnijrn^S m ) ^ m _i,i m _i,5 m _i(j m _i) \^m— 1) jm— 1) £>m— l) ■ (2) 

These equations use the requirement that D be a local operator in that its values depend 
only on the state of the qubit at the head location. It is independent of the states of qubits 
at other lattice sites. In this case one has D\l,j,S) = Di^ t s\l,j, S) = Di t g(j)\l,j,S). The 
connection between the eigenvalues of D and the values of 7; iS is then given directly by 

D l,S(j) = ll,S(j)- 

If T and are distinct path generating on B then the direct sum decomposition of T 
means that the Feynman Hamiltonian [40 of Eq. [I] is also decomposable as H = J2i HiPi 



where for each i, 

H l = K(2-T l -Tl). (3) 

Since T« is a weighted shift on the path of states in Bi, Hi is a tight binding Hamiltonian 
on Bi and thus in the subspace PiH. 

As is noted elsewhere [^] Hi can also be written in the form 

Hi = K(2 -Ui- Uj) + Vi (4) 

where K(2 — Ui — Uj) is, formally, the kinetic energy associated with the evolution of the 
computation along path i. For the physical model used here it is also the kinetic energy of 
head motion on the qubit lattice. (The kinetic energy depends on the magnitude but not 
the direction of the head momentum on the lattice.) 

The potential V- = K(U { -Ti + Uj - T}) is a nearest neighbor off-diagonal potential with 
matrix elements on internal path states given by 
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(mf, i\Vi\m, i) = K[(l - A,m)<W,m+l + (1 ~ A,m-l)5m',m-l]- (5) 

Since |m,z) = \l m ,jm, 3 m ) and D is local, D i>m = 7z m ,s m (j m )- 

In Schrodinger evolution under the action of H, Eq. [1], the choice of which Hi is active 
is determined by the initial state \P(0). If \l/(0) is a superposition of of computation basis 
states in just one path i, Then Hi is the only active component. If ^(0) is a superposition 
of computation basis states in different paths, then all Hi are active for just those paths i 
in which \1/ (0) has a nonzero component. Note that in this case the states in the different 
paths evolve coherently. 

The simplest nontrivial choice for D (other than D = 1) is that the D has two eigenvalues 
1,7 where 7 is a positive real number between and 1. This will be implemented here by 
limiting consideration to generalized quantum Turing machines (GQTMs) that have a fixed 
potential associated with the read-1 elementary steps only. All other qubit read steps are 
potential free. Read-1 steps are defined to be all steps \m,i) — > \m + in which an 
elementary step term 7]^ is active (i.e. the qubit at the location of the head is in state |1)). 
This is accounted for by setting 

r = £[(£Wi,-)+7Wi,i] (6) 

I S^l 



IV. THE COUNTING GQTM 

A. The Step Operator 

The counting GQTM serves as an interesting example of the foregoing. Based on Eq. |5| 
the step operator for this GQTM is a sum over seven elements or terms (each term includes 
the j-sum): 

00 

T = £ {QoPojuPj + wQoP^uPj + Q^ojuPj 

j=-oo 

+ wQxP 2 jU^Pj + lQ2V x jPljU ] Pj + n A () 1 r ri l\ ]j iil' J 

+ wQ 2 P 23 uP 3 ) (7) 

The projection operators Qi, P s j Pj are for the head in state |Z), the site j qubit in state \s), 
and the head at site j. w is a shift mod 3 on the three head states (wQ m = Q m+ iw mod 3) 
and u shifts the head along the lattice by one site (uPj = Pj + iu). The need for markers is 
accounted for here by choosing the qubits in the lattice to be ternary with states |0), |1), |2). 
\2) is used as a marker and |0), |1) are used for binary strings. The qubit transformation 
operator v X j = cr X j(P j + Py) + P 2 j exchanges the states |0), |1) and does nothing to the 
state 1 2) for the site j qubit. 

The adjoint is defined from T in the usual way noting that the operators for the 
head states, head position states, and qubit states commute with one another. Note that 
v lj — v xj- This gives 
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j=-oo 

+ Q 1 w ] P 2j P j u + ^Q2P\jV X jPjU + Q 2 wP jV xj Pj 



u 



t 



+ Q 2 w ] P 2j P j v)) (8) 

The terms in T have been chosen so that, besides carrying out the desired operations, T 
is distinct path generating in the computation basis. As defined T and satisfy a sufficient 
condition for this [^7| , namely, that T^T = J2t=i an d TT^ = J2t=i ^ where t and 



denote terms of T and in Eqs. [7| and |[ 

For this example all steps are potential free except those in which the 5th term is active. 
For these a potential, given by Eq. [5], is present as D^i = 7 for I = 1. The distribution of 
potential and potential free steps along a path depends on which of the 7 terms is active for 
a path state. The potential width at any path location is given by the number of successive 
iterations for which the 5th term is active. 

Comparison with Eq. || shows that for this GQTM (T = UD) U is a also sum of 7 terms 
with C/2,1 the fifth term of Eq. |7| with 7 excluded. The other 6 terms of U in the I, s sum 
are given directly by Eq. ^. 

To apply this to a specific example, consider the initial state, shown in Figure 1, with the 
head in state |0) in a wave packet localized to the left of the origin. All qubits are in state 
|0) except those at sites 0,n + 1 which are in state |2). The initial head and lattice qubit 
state are the path labels for this example. Iteration of T on this state generates in turn 
all the integers as binary strings of length < n. When the space between the two markers 
is completely filled with Is, corresponding to the integer 2 n — 1, the last pass of the head 
changes all Is to Os. The head in state |1) then moves to the right away from the marker 
region as the enumeration is completed. A more detailed description of this process, based 
on iteration of T, is given elsewhere ||22||. 



B. Potential Distributions 

The path potential distribution for the counting GQTM is closely related to the distri- 
bution of read-1 steps in all the steps (read-1, read- 0, and read-2) obtained by iteration of 
T. This distribution can be represented by a function TZ from the nonnegative integers to 
{0, 1} where TZ(j) = 1[0] if the jth step (or T iteration) is a read-1 [read-0 or read-2] step. 
The initial state or origin is assumed to be with the head in state |1) at position n+ 1. The 
location of the origin is not important provided all read-1 steps are included. 

The read-1 steps occur during a stepwise search, starting from the units place, for the 
first after a string of Is. Based on this TZ would be expected to be related to the function 
R from the nonnegative integers to the nonnegative integers such that, in the binary string 
for j, R(j) gives the number of Is occurring before the first 0. Because of this, the properties 
of R can be used to obtain a useful expression for TZ. 

The range sequence of R, which is 0,1,0,2,0,1,0,3,-- - , is known as the hierarchical 
sequence It can be expressed recursively [^2| by the pair of recursion equations; Rn = 



Sn^iN, Sn = -RatStv-i for iV = 1, 2, • • • with So = 0. R^ is the first 2 N elements of R. R is 
a generalized substitution sequence under the substitution rule n On+1 for n = 0, 1, 2, • • •. 



7 



The generalization denotes the fact that, contrary to the literature definitions fl4,42|)43|, the 
alphabet is infinite. 

The recursion relations for R also hold for 1Z if N is replaced with N_ where N_ denotes 
the string of numbers Ol^O^" 1-1 . Here 1^ and N+1 denote strings of N Is and N + 1 Os. 
That is, for N = 1, 2, • • • 

TZn = Sn-iN_; Sn = TZ-nSn-i- (9) 

where S = 0. 

This equation gives the distribution of read-1 steps in all steps associated with enumer- 
ating the first 2 N integers. An example for iV = 6 is given elsewhere [22|. The explicit 
locations of bands of read-1 steps can be obtained by noting that the number of elements in 
Sjy and Tl^ is given respectively by 2 N+3 — 2(N + 3) and 2 N+2 — 2. The number of read-1 
bands is equal to 2 7V_1 . 

The relation of the read-1 step distribution to the potential distribution in the tight 
binding Hamiltonian, Eq. [|, for the path corresponding to enumerating the first 2 N integers, 
can be obtained from TZ^f. Extend TZn to a two way infinite sequence by adding 0s to both 
ends. Define a function Di on the spectrum of D by D^ m = Di(m). The function D{ is 
related to the extended 7Z n by 

A = 1 + (7 - l)K N . (10) 

. From this the matrix elements of Vi , Eq. |5|, are given by 

(m'\Vi\m) = K(l - j)[fi N (m)5 m > tm+1 + TZ N (m - l)5 m ', m -i]- (H) 

Eq. |TT] shows that the off-diagonal potential barrier associated with a band of n Is at 
path sites a, • • ■ , a + n — 1 with 0s on both sides (i.e. lZjsr{j) = 1 if a<j<a + n — 1 and 
H N (a — 1) = TZjy(a + n) — 0) extends from sites atoo + n. The potential barrier has a core 
of n — 1 sites of height 2K(1 — 7) flanked by single site potentials of height K (1 — 7) at sites 
a and a + n. This is the case whether these sites correspond to initial or final states in the 
matrix elements of Eq. [11]. Note that a potential is associated with the site a + n for which 
TZn has value 0. If n = 1 the barrier is two sites wide and has a height of K{1 — 7) as no 



core is present. The potential distribution given by Eq. [XT] for iV = 6 is shown in Figure 2. 

Figure 2 is a good illustration of the general property of the distribution, that barriers 
that are two sites wide (n = 1) are the most common as 1/2 the barriers are of this type. 
One half the remaining barriers are 3 sites wide n = 2, and 1/2 those remaining are 4 sites 
wide (n = 3), etc.. 

There are many other possible initial states and paths with associated tight binding 
Hamiltonians besides that given in Figure 1. Discussion of these, and the collection of paths 



into equivalence classes, is discussed elsewhere [22|. 



C. Transfer Matrix 



Following Erdos and Herndon [18], solutions for tight binding Hamiltonians, such as that 
riven by Eq. £|, can be written as — + + ^ ni where and ^ ni are the states 
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coresponding to the path boundary potential free regions and is the state in the path 
region containing the potentials. For example if the path regions for j < and j > n are 
potential free (integers label path states) then 



-i 



tf 7 = Y, (Ae ikj + Be- ikj )\j) 



j=-oo 
oo 



j=n+l 



(12) 



Here j is the path state label and k is the momentum in the potential free regions. 

The form of depends on the form of the potential. For a single barrier corresponding 
to a band of N read-1 sites extending from site a to a + N — 1, 



a+N-l 



*// = {Cie iha + D ie - iha )\a) + £ (C 2 e ilj + D 2 e~ ilj )\ 3 ) 

j=a+l 

+ (C 3 e ih{a+N) + D 3 e- ih{a+N) )\a + N). 



(13) 



If the energy is below the barrier height of 2K(1 — 7), which is the case of interest here, 
then I is replaced by il. If the energy is below K(l — 7) then the momentum h is replaced 
by ih. If several potentials are present is a sum of terms of the form of Eq. [13|. 

The four complex coefficients A,B,F,G are completely determined by boundary con- 
ditions for the case under study and the properties of the transfer matrix. For example if 
initial states are required to be localized in region I (Figure 1), then G = as there is no 
incoming head motion in region III. 



The transfer matrix Z relates the coefficients F, G to A, B according to [18 



Z is a unimodular (determinant=l) 2x2 matrix which also satisfies [|l| 



Z\2 



^22 



J 21- 



(14) 



For a string of potentials the transfer matrix connecting the state on the right to the 
state on the left is the product of the transfer matrices for each of the potentials |18|-|2"0|. 
If the potential distribution is given by a recursion relation, such as that of Eq. ||, then 
the transfer matrix Z^ for the potential distribution corresponding to TZ n is given by the 
recursion relations 



Zn — WnXn-i 

X N = X n _ 1 WnX n _i 



(15) 



with X |n = e and X |i2 = 0. Eq. |T3| holds for X . is the transfer matrix for the 
potential corresponding to the sequence N_ = 01 N N+l . As noted earlier this corresponds 
to a potential barrier with a core width of N — 1 sites flanked on each side by a one-site 
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potential half the core height. Potential free regions of 1 site and iV + 1 sites surround the 



barrier. Note that the order of matrix multiplication in Eq. [15] is the inverse of the order in 
which terms appear in the sequence of Eq. |9|. 

This use of recursion relations has the advantage that in order to obtain the elements 
of the matrix Zn, only polynomially many matrix multiplications are needed. If Zn is 
obtained from the matrices associated with each potential barrier in the distribution, then 
exponentially many matrix multiplications are required. 

From Eq. [15] one sees that the only matrices needed in explicit form (other than X ) are 
the W m for m < N. These can be obtained explicitly (see also (|TJ|) by setting (E — Hi)^/ = 



where Hi is given by Eqs. [| and 11, and the state \1/ is given by Eqs. O and [L3[ A set of 



linear equations is obtained by setting the coefficients of each state \j, i) equal to 0. that is 

y EV,i\E-H i \j',i){j'\*) = Q. (16) 

3' 

Because Hi contains nearest neighbor interactions only, just 3 terms in the j' sum, j' = j, j±l 
contribute. 

The linear equations so obtained can be used to derive the matrix elements of W m for 
m = 01 m m+1 . A brief summary of the derivation is given in the Appendix. For energies 
below the potential core height (i.e. less than 2K(1 — 7)), they are 



gjfc(m+2) 

W m \ u = : — [e 2ifc sinh/m - 27e ifc sinh/(m - 1) + 7 2 sinh Urn - 2)1 

2ry sin k sinh I 



W m \i2 = : —[sinh Im — 27 sinh Urn — 1) cos k + 7 2 sinh Um — 2)1 (17) 

2ry sin k sinh I 

where the momenta l,k are related by 2K(1 — cos A;) = E = 2K(1 — 7 cosh I). For the 
flank potentials or those for m = 1, E = K{2 — (•y 2 + 27cos2/i + l) 1 ^ 2 ). However, the 



momentum h, Eq. [13|, does not appear in the above. These equations can be shown by 
explicit calculation to hold for m — 1 and m = 2 as well as for larger m. One also has 
W m \ u = W^\ 22 , W m \ 2 i = W^\i2] and W m is unimodular. 



V. ENERGY BANDS AND GAPS 



As is well known p5jJT9|j2l]| ] the distributions of energy bands and gaps for the infinite ID 



periodic system, with a unit cell satisfying the problem being considered, can be obtained 
from a plot of TrZ^ as a function of the momentum k. Regions for which \TrZ^\ < 2 
correspond to energy bands; regions with \TrZ N \ > 2 correspond to energy gaps. As the 
main interest in the paper is in bound states, k is restricted to be less than arccos7. The 
infinite periodic system is referred to as the associated periodic system. 

Test calculations show that except for values of 7 very close to 1, the values of TrZ^ 
fluctuate extremely rapidly between astronomically large positive and negative values (of the 
order of ±10 where M is very dependent on 7 and N. For values of 7 appreciably different 
from 1 and values of N of order 10, M has values greater than 200. It is thus necessary 
to determine values of the two free parameters K, 7 that are reasonable physically. To do 
this it is necessary to assign a spacing A to the lattice qubits. Then the dimensionless 
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momenta /, k are related to the usual momenta I', k' by k = k'A, I = I' A. The form of K 
is obtained from the fact that if T = Y, the bilateral shift on the lattice, then H, Eq. [j], is 
the symmetrized form of the second derivative on the lattice. This gives 

h 2 c 2 

K = ^~, (18) 
2mA 2 v ; 

where m is the mass of the head in energy units (mc 2 ). 
Since V = 2K(1 - 7), with 7 = 1 - 5, 



r 2mVA 2 

' — »ps- (19) 



For electron systems, m is of the order of the electron mass, V is a few electron volts, and 
A is measured in Angstroms. Taking m equal to 2 electron masses, A = 1A, and V = 2ev 
gives S ~ 0.001 or 7 ~ 0.999 

The method used for calculation of TrZ^ consists of using the recursion relations, Eq. 
1~5], to obtain and then taking the trace. The well known method of using trace maps 



42|,PJT2]|, is not used because it is more efficient to obtain directly the matrix and 
then take the trace rather than by use of the trace map relations. In particular, obtaining 
Zn directly from Eq. ITR using Eq. [14], and then taking the trace requires a number of 
scalar multiplications proportional to N. Use of the trace map requires a number of scalar 
multiplications proportional to N 2 . 

Energy bands and gaps were calculated as a function of the dimensionless momentum 
k for several different values of 7 and N. The results for 7 = 0.999 are given in Figure 3 
for iV = 10, 15, 18 as three bands of vertical lines connecting horizontal line segments or 
points. The upper or lower horizontal line segments or points represent respectively energy 
bands or energy gaps. Points represent bands or gaps whose width is of the order of the 
values of the increment or spacing AA; used to make the calculations. 

The range of momentum shown is from k = 0.0223 to k — arccos7. The lower limit is 
chosen because, to the accuracy of the calculations, there are no energy bands in the region 
k < 0.0223 in which TrZjq increases rapidly to very large values. Each band shows the 
result of calculations using about 2800 equally spaced grid points which gives Ak « 7/56OO. 

The lower band for iV = 10 shows that most of the momentum region 0.0223 < k < 
arccos7 = 0.0447 is covered by energy bands separated in most cases by very narrow gaps. 
The widest gap is at k ~ 0.033 with smaller gaps at k ~ 0.025, 0.029, 0.038, and 0.043. 

The middle band for N = 15 shows that the wide energy bands for N = 10 are 
split into many smaller bands with intervening gaps. Gaps which were already present 
at iV = 10 continue to be present. This splitting continues for increasing N as is 
shown by the upper band for N = 18. The width and positions of the larger gaps at 
k ~ 0.043, 0.038, 0.033, 0.029, 0.025 are maintained. The bands also show that as N 
increases the number of bands increases greatly. However the band widths decrease even 
faster so that the fraction of the momentum region associated with energy bands, although 
sitll large, decreases as N increases. 

The distribution of bands or band density is shown by the spacing or density of the 
vertical lines. For example for iV = 18 the density is high in the dark regions of the band, 
such as around k w 0.024 and 0.036. It is low in the regions of the wide gaps. 
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Calculations were not made for values of iV > 18 because the range and fluctuation of 
values of TrZ^ becomes so extreme that the accuracy of calculations that can be made 
with a reasonable amount of time and effort becomes questionable. Even at N = 18 high 
resolution calculations (to be discussed later) over small momentum regions are needed 
to resolve individual bands. Recall that the results shown give the energy band and gap 
structure for an infinite periodic structure with a unit cell of 2 N ~ l potentials with a width 
and location distribution given by Eq. [| Thus for N = 18, about 128, 000 potentials are 
present in the unit cell. 

In all figures except Figure 1, the values of the momentum used in the abcissa are 
larger than the actual values by a factor of 1000. They are also dimensionless. This has the 
advantage that the bands in Figure 3 are valid for any values of m, V, A that give 6 = 0.001 
(Eq. [19]) or 7 = 0.999. One set is the values chosen for the electron range. Another set 
in the nuclear range (probably impossible to realize) that gives 5 ~ 0.001 is m equal to 
the proton mass, V = lMev and A = 3F. (IF = 10~ 13 cm). For the values chosen in the 
electron range the momentum range shown in the figure corresponds to an energy range of 
from 1 to lev . 

It is useful to explore a wider variation of parameters, such as larger values of m, A, 
or V which give different values of 7. Choosing values which increase values of 5 gives one 
more latitude in a choice of physical model parameters. For this reason calculations were 
also made for 5 larger by a factor of 10. The results are shown in Figure 4 for 7 = 0.99 for 
the same values of N used in Figure 3. The momentum range is cut off at k = 0.069 as 
there do not appear to be any bands below the cutoff. 

Compared to the results for 7 = 0.999, Figure 4 shows a great increase in the fraction 
of the momentum region occupied by energy gaps. There are large gaps around values of 
k = 0.084 and 0.10 and especially 0.12. There are many more energy bands present than 
are shown in Figure 3 for the same value of N. But the widths are so much smaller that the 
fraction of the region alloted to energy bands is much less. This is especially pronounced in 
the bands for N = 15 and N — 18 where energy bands show as points only. For these values 
of N, the maximum width of any energy band is only about 0.0002. The calculations for 
N = 15 and N = 18 should only be taken as giving a rough indication of the bands present. 
Use of finer grids in the calculations would be expected to show more bands of very small 
widths in the high density band regions than are shown in the figure. The large gaps would 
be expected to remain. 

Increasing the value of S by another factor of 10 to give 7 = 0.9 gives the results shown in 
Figure 5. They were obtained using a grid with over 7000 equally spaced points with a low 
momentum cutoff at 0.166. The values of N used are smaller than those in the preceding 
figures because the fluctuations in the values of TrZ^ are so large. 

The calculations show that very few energy bands are present and except for N = 6 the 
widths are pointlike. For N = Q the band grouping appears to be hierararchical in that pairs 
of doublets are grouped into pairs of quartets which are in turn grouped into octets. The 
appearance of the band at k ~ 0.196 as a singlet instead of a doublet may be a consequence 
of the finite resolution of the calculations. The gaps between the groups increase as the 
group size increases. There are indications that this structure may be preserved for N = 8 
and possibly for N = 10. However, each band splits up into many smaller bands as N 
increases. More detailed calculations with extremely fine grids would be needed to check 
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this. 

The reason for the low momentum cutoffs in the appearance of bands is not clear. They 
do not appear to be associated with the height of the narrowest barriers (1/2 of the total),i.e. 
those with height K(l — 7). For these barriers the corresponding values of the momentum 
at which the energy equals the barrier height (cos A; = (1 +j)/2) are 0.0316, 0.100, 0.317 
for 7 = 0.999, 0.99, 0.9 respectively. These values are appreciably higher than the values 
at which bands first appear. 



VI. LANDAUER RESISTANCE 



Interest in the Landauer Resistance [17,18], defined as the ratio of the reflection prob- 



ability to the transmission probability through a finite sequence of potentials, is based on 
the fact that for ID electronic systems, it is a measure of the transmission of the electron 
through the system. 

From the standpoint of this paper the Landauer Resistance (LR) is of interest because the 
reciprocal, 1/(1 + Li?), is the probability of transmission through the potential sequence. For 
any path of states for a GQTM, the reciprocal is a measure of the probability of completion 
of the quantum computation as it refers to the transmission of a wave packet from one end of 
the path to the other through potentials associated with the path states. For the counting 
GQTM it is a measure of the probability of completion of the counting with potentials 
associated with read-1 steps only. 



Landauer |23| has emphasized that reflection from background potentials associated with 
each step of a quantum computation degrade the computation due to multiple reflections 
from and transmissions through the potentials. In particular for energies below the potential 
height the transmission would be expected to decay exponentially as the packet must tunnel 
through each of the potentials, 

There are indications that this may not be true for potential sequences which are not 



random. For example it has been shown [ IS , 2C ] that for Thue- Morse, Fibonnaci, and periodic 



sequences in the Kronig-Penney model, the Landauer Resistance fluctuates rapidly with 
energy between very high values >> 1 and very low values << 1. In addition high values 
are associated with energy gaps in the infinite associated periodic system; low values are 
associated with energy bands. 

For these reasons the Landauer Resistance (LR) has been calculated for the counting 
GQTM for various values of N and 7 as a function of the dimensionless momentum k. The 
calculations are done using the result that the LR is given by [JIT 



LR=\Z N \ 12 \ 2 (20) 

where is given by Eq. [15]. The results are shown in figures that show the functional 
dependence of the log (base 10) of the Landauer Resistance (log LR) on the momentum. 
The momentum units are dimensionless. 

The results for 7 = 0.999 and iV = 10 are shown in Figure 6. The enery band and 
gap structure for the associated periodic system is reproduced at the bottom of the figure. 
The calculations show a rapidly fluctuating log LR with large maxima at momentum values 
corresponding to energy gaps. The heights of these large maxima appear to be correlated 
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with associated gap widths; wide gaps have greater maxima than narrower ones. Also the 
general trend of the values such as the height of the maxima, is towards lower values as the 
momentum increases. Both these features were noted by Roy and Khan [|T9| , |20fl . However 
for the smaller maxima, which occur in regions of energy bands of fairly uniform width 
separated by narrow gaps, the opposite appears to be the case in that the maxima appear 
to be associated with energy bands. In addition log LR has sharp spikes to very low values 
at the regions of these narrow gaps. 

An interesting feature of the results is that for most values of the momentum log LR 
is less than 0. A value of corresponds to a 50% probability of transmission through the 
potentials. Thus for this case the probability of completion of the quantum computation 
is > 50% for most momenta. For selected values the probability is close to 100%. These 
low values are very likely correlated with the result that most momenta are in energy bands 
with only a small fraction of the momentum region taken up with gaps. 

Figure 7 shows a higher resolution plot of log LR for a portion of the momentum region. 
It shows very clearly the association of the sharp minima to band edges. A study of the 
actual computation points for the figure shows that the minima are associated with the 
upper edge of a band just below a gap and not with the lower edge of a band just above a 
gap. Also the depth of the minima appear to depend on the width of the associated gap in 
that in general the depths are much less for wider gaps than for narrower ones. For example 
the minimum with the relatively wide gap at k — 0.0355 is higher than the sharp spike 
minima on either side with relatively narrow gaps. The minimum associated with the even 
wider gap at k = 0.038 has almost disappeared. 

The results of a calculation of log LR for the same value of 7 but for a larger value of 
N = 18, are shown in Figure 8. Just as the band spectrum is much more finely divided 
log LR fluctuates much more rapidly and violently. For the purposes of clarity the plot is 
cut off at values of log LR > 20 (LR > 10 20 ). The main reason is that the interest here is 
in the behaviour of log LR around and not how high the values can go. 

The results show that there are a great number of values of the momentum for which 
log LR is or less. However there are also a large number of fluctuations to much higher 
values. Because of the finely divided nature of the spectrum and extreme violence of the 
fluctuations, the figure should be taken only as an approximate indication of the dependence 
of log LR on the momentum. 

A higher resolution calculation over the range 0.032 < k < 0.033, Figure 9, shows very 
clearly the preponderance of values at or below 0. as well as the rapid fluctuation to large 
maxima. The heights of the maxima are clearly correlated with the widths of the associated 
gaps. However, the individual small fluctuations are barely resolved. The heights of the 
peaks also show a regularity that can be correlated with the widths of the potentials in the 
distribution. The regularity appears to be correlated with the hierarchical structure of the 
sequence given by Eq. |9|. 

A very high resolution calcuation over a very small momentum region of 0.03225 < k < 
0.03228 for 7 = 0.999 iV = 18 is shown in Figure 10. The figure shows very clearly the 
correlation between the decrease in the depths of the minima and the increase in the widths 
of the associated gaps in the band spectra. For example the gap at k = 0.032261 is sufficiently 
wide that the associated minimum in log LR has just about disappeared with only a slight 
dip present. Shallow minima are associated with the narrower gaps at k = 0.032253 and 
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k = 0.032269. The figure shows also that band widths are much more uniform than gap 
widths which vary widely. This point was noted before for smaller values of N. 

These figures show that there are fairly wide regions of the momentum where, except 
for a few sharp spikes, log LR < 0. It is noteworthy that this corresponds to transmission 
through more than 128, 000 potential barriers with the energy below the barrier height for 
1/2 these potentials for 0.0316 < k < arccos7. For k < 0.0316 the energy is below the 
barrier height for all the potentials. It is possible that the high transmission is due to 
the quasiperiodicity of the potential distribution; however more work is needed to test this 
possibility. 

In order to show the dependence of log LR on 7 for fixed N, calculations were made 
for iV = 10 and 7 = 0.99. this value is larger by a factor of 10 than that used so far. It 
corresponds to the electron regime with larger parameters in Eq. |1£] to give 5 = 0.01. The 
results are shown in Figure 11. 

The figure shows very clearly the correspondence between maxima in log LR and gaps in 
the spectrum. Large gaps are associated with high peaks and small gaps with small peaks. 
Also the correlation between the peak heights and the hierarchical structure of TZ, Eq. |9|, 
appears to be present. In regions where the band and gap widths are small, log LR fluctuates 
between values of about 1 and —2 or less. On average the values of log LR, even in small 
gap regions, are higher than for 7 = 0.999. This is to be expected as smaller values of 7 
correspond to larger values of the potential. On average the transmissivity is appreciable as 
there are extended regions where the values of log LR are around 0. 

Calculations for S larger by another factor of 10 to 5 = 0.1 (7 = 0.9) are shown in Figure 
12 for N = 8. The extreme violence of fluctuations in the values of log LR, even for the 
relatively low value of iV = 8 is evident. This is to be expected since higher values of 5 
corresponds to higher values of the potential (with m and A fixed). Also since most of the 
space is occupied by gaps, log Li? has very high values and approaches only in the regions 
of narrow bands. 

A high resolution calculation for 7 = 0.9 and N = 8 over a small portion of the mo- 
mentum region 0.40 < k < 0.415 is shown in Figure 13. The figure shows very clearly the 
association of minima in the LR with upper edges of energy bands and the loss of minima 
with increasing gap width. The gap including the value of k — 0.408 is so wide that the 
associated minimum in log LR is completely gone. Also the values of the LR in the small 
momentum region shown in the figure are between 10 and 100 on average with values at or 
below 1 (log Li? < 0) occuring only rarely. 

VII. DISCUSSION 

To summarize, the main characteristics of the band spectra shown here include the 
relative uniformity of band widths compared to gap widths over the bound state momentum 
regions. This is especially apparent in Figure 3. This figure and Figure 4 suggest that 
fluctuations in bands widths do appear to increase somewhat as 5 increases. 

The figures also show that as 5 increases the gap widths and the fraction of momentum 
range occupied with gaps increases greatly. For 5 = 0.001 most of the k region is occupied 
by bands. For 5 = 0.01 a large fraction is taken up by gaps and for 5 = 0.1 almost all the 
region is gaps with a few very narrow bands present. 
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Figures 3-5 show that for fixed values of 7 as N is increased, energy bands are split into 
increasing numbers of bands with decreasing widths. Gaps that exist for a given value of N 
appear to be preserved as N is increased above the given value. This suggests but does not 
prove that in the limit of N = 00, for each 7 the spectrum of the tight binding Hamiltonain 
is singular continuous and is a Cantor set. However, the proofs in the literature do 
not hold here as the generalized substitution sequence of Eq. || is not primitive. 

The values of log LR also fluctuate widely over the momentum range with the fluctuations 
very dependent on the values of 5 and N. Sharp minima are located at or near band-gap 
boundaries with more minima located at or near gap-above-band edges than at or near band- 
above-gap edges. The reason for the location of minima at or near band gap boundaries is 
not clear. It is worth noting that, based on the properties of Zjy, Eq. LR = \Zn\i,i\ 2 — 1. 
At the boundaries \TrZ N /2| = 1. This gives 

LR=\Q(Z N \ ltl )\* (21) 

at the boundaries. 5s (x) denotes the imaginary part of x. 

Another aspect of the minima is that none are associated with boundaries of sufficiently 
wide gaps. The figures show that in almost all cases the depth of the minimum in log LR 
is very dependent on the width of the associated gap. If the gap is narrow, the minimum 
is deep. As the gap becomes wider the associated minimum becomes shallower until it 
disappears for sufficiently wide gaps. An example of this where the minima has almost 
disappeared is the gap just below k = 0.03227 in Fig 10. Other examples can be found in 
the other figures. 

One consequence is that for wide gaps the maxima in the LR appear to be centered over 
the gaps. For narrow gaps there is an associated minimum in the LR . For these cases the 
maxima in the LR appear to be centered over the bands and not the gaps. 

These results differ from those obtained by Roy and Kahn [frU|,pO| for the Thue-Morse 



lattice in that they found all maxima in the LR centered on the gaps irrespective of the gap 
widths. The reason for this difference is not clear. The results obtained here do agree with 
those of Roy and Khan in that the larger LR maxima are associated with wide gaps, and 
the smaller maxima are associated with narrower gaps. As noted, here the smaller maxima 
tend to be centered over the bands and not the gaps. 

As is well known the Landauer Resistance is a measure of the transmission through a 
sequence of potentials. For the case at hand it is a measure of the completion probability for 
the counting generalized quantum Turing machine to enumerate the nonnegative integers 
up to 2^ — 1. However the strong dependence of the LR on k means that care must be taken 
in its use. For example Figure 1 shows an initial state represented as a wave packet over 
different head lattice positions with a fixed lattice qubit state. Transformed to momentum 
space the initial state is in essence a wave packet over the head momentum k with the same 
fixed head internal state and lattice qubit state. That is \& = JZ n Ck\l, k, S)dk where |/, k, S) 
is a state for the head in state I and momentum k and the qubit lattice in state \S). 

Linearity of the Schrodinger Equation gives the result that, for each component state 
Ck\l, k, S), the transmission and reflection coefficients are given respectively by 

p _ °k 

2,2 
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B = (22) 

Zn\22 

Here Eqs. [12] and [H| were used along with the properties of Zn, Eq. [14]. The coefficient 
G = as there is no incoming component in region III. As the results in Section |V| show, the 
coefficients F and B are very dependent on k. Of course the LR, given by LR = \B\ 2 /\F 2 \, 
is independent of the normalization to shown in Eq. ^2[ 

More general linear superpositions in the initial state are possible. For example the most 
general initial state for the counting GQTM is given by ^ = J2i,s J-n ^> ^? Each 
component in the S and I sums corresponds to a different input state and a different quantum 
computation. This ability of quantum computers to carry out simultaneously computations 
on different input states i.e. parallel quantum computation was first described by Deutsch 
2f . 



From the perspective of this paper an interesting aspect of the above is that Schrodinger 
evolution with one Hamiltonian, Eq. |l] for the system is equivalent to a linear superposition 
of different Schrodinger evolutions each with a different tight binding Hamiltonian. This 
follows from the fact that in general for each component in the I, S sum, the potential 
distribution on the corresponding state path depends on \S) and on |Z). This is an example 
of a linear superposition of the action of component Hamiltonians in the decomposition of 
the Feynman Hamiltonian, Eq, |l| into a sum over tight binding Hamiltonians, Eq. [3]. 
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APPENDIX 

As noted in the text around Eq. |TJ| solution of (E — H)^/ = gives a set of linear equations, 
one for each path state In regions of constant potential or that are potential free, 7 = 1, 
the linear equations have the same form: 

(jXY)(E - 2K) + Kj(j - 1XY) + K*y(j + IXY) = 0. (23) 

Here (jXY) is a shorthand notation for Xe iqj + Ye~ iqj where A, Y = A, B; C*, D { \ or F, G 



respectively of Eqs. 12 and 13 and q is the appropriate momentum. 



Solution of this equation gives E = 2K(1 — 7cosg). For a band of m read-1 sites (from 



a to a + m — 1) surrounded by 0s on both sides, use of Eqs. [L2|, |13, and [16] gives 6 equations 
different from Eq. for the sites a — 1, a, a + 1, a + m — 1, a + m, a + m + 1. In order of 
these sites, they are 

(aAB) - {aC x D x ) = 
(aCi£>i)(£ - 2K) + K((a - 1AB) + A 7 (a + IC 2 D 2 ) = 
(a + IC 2 D 2 ){E - 2K) + K 1 (aC 1 D 1 ) + A 7 (a + 2C 2 D 2 ) = 
(a + m - 1C 2 D 2 )(E - 2K) + A 7 (a + m- 2C 2 D 2 ) + A 7 (a + mC 3 D 3 ) = 
(a + mC 3 D 3 )(E - 2K) + A 7 (a + m- 1C 2 D 2 ) + K(a + m + 1FG) = 

(a + mFG) - (a + mC 3 D 3 ) = 0. 
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Here Eq. ^ has been used to simplify the first and last equations to remove (a — 1AB){E — 
2K) and (a + m + 1FG)(E - IK). 

The first equation is used to remove (aCiDi) from the next two equations. Combination 
of the resulting two equations with Eq. ^ gives a 2 x 2 matrix connecting (a + IC2-D2) to 
(aAB). In a similar way the second 3 equations are used to give a 2 x 2 matrix connecting 
(a + mFG) to (a + m — IC2D2). The product (taken in the correct order) of these two 
matrices and a matrix connecting (a + m — IC2-D2) to (a + IC2-D2) gives the final result. 
The matrix for W m is obtained by setting a = 1 and translating the result by another m + 1 
potential free sites. 



18 



REFERENCES 



[1] D. Shechtman, I. Blech, D. Gratias, and J. Cahn, Phys. Rev. Letters 53 1951 (1984). 
[2] T. Fujiwara and T. Ogawa (Eds) Quasicrystals, Proceedings of the 112th Taniguchi 

Symposium, Shima, Mie Prefecture, Japan Nov 14-19 1989, Springer Series in Solid 

State Sciences 93, Springer Verlag New York. 
[3] D. DiVincenzo and P. Steinhardt Eds., Quasicrystals the State of the Art, in Directions 

in Condensed Matter Physics- Vol. 11, World Scientific Singapore 1991. 
[4] A. Bovier and J-M. Ghez, J. Phys. A: Math. Gen. 28 2313 (1995). 
[5] Let A be a finite alphabet, A* the set of all finite words of A, and A N the set of all 

one way infinite words of A. A substitution rule replaces each letter a of A with a word 

from A*. For each a a sequence of words in A* can be generated by starting with a and 

iterating the substitution rule. The limit sequence in A N is a substitution sequence if it 

is invariant under the substitution. 
[6] J.M. Luck, Phys. Rev. B 39 5834, (1989) 

[7] J. Bellissard, A. Bovier, and J-M. Ghez, Commun. Math. Phys. 135 379 (1991) 
[8] A. SiitS, Jour Stat. Phys. 56 525 (1989). 

[9] A. Bovier and J-M. Ghez, Commun. Math. Phys. 158 45 (1993); 166 431 (1994). 
[10] M. Hornquist and M. Johansson, J. Phys. A Math. Gen. 28 479 (1995). 
[11] C. S. Ryu, In-mook Kim, G. Y. Oh, and M. H. Lee, Phys. Rev. B 49 14991 (1994). 
[12] Y. Avishai, D. Berend, and D. Glaubman, Phys. Rev. Lett. 72 1842 (1994). 
[13] Y. Avishai and D. Berend, Phys. Rev B-II 43 6873 (1994). 
[14] Y. Avishai and D. Berend, Phys. Rev B-II 45 2717 (1992). 

[15] V. M. Gasparian and A. Gh. Khachatrian, Solid State Commun. 85 1061 (1993). 

[16] P. Erdos and R. C. Herndon, Solid State Commun. 98 495 (1996). 

[17] R. Landauer, Phil. Mag. 21 863, (1970). 

[18] P. Erdos and R. C. Herndon, Adv. Phys, 31 65, (1982) 

[19] C. Roy and A. Khan, Phys. Rev. B, 49 14979 (1994): 

[20] C. Roy and A. Khan, Solid State Commun. 92 241 (1994). 

[21] C. Roy and C. Basu, Phys. Stat. Sol. (b) 177 315 (1993). 

[22] P. Benioff, To Appear in Phys Rev. Letters 

[23] R. Landauer, Philos. Trans. Roy. Soc. London A353 367 91995); Physics Letters A 217 

188 (1996); Physics Today 44 No. 5, 23 (1991). 
[24] A. Peres, Phys. Rev. A32 3266 (1985) 
[25] W. H. Zurek, Phys. Rev. Lett. 53 391 (1984) 

[26] P. Benioff, Jour. Stat. Phys. 22 563 (1980); Ann. NY Acad. Sci. 480 475 (1986). 
[27] P. Benioff, Phys. Rev A 54 1106 (1996). 

[28] D. Deutsch, Proc. Roy. Soc. London Ser. A 400 97 (1985); 425 73 (1989). 
[29] S. Lloyd, Sci. American 273 No.4, 140, (1995); Journal of Modern Optics, 41 2503 
(1994). 

[30] D. P. DiVincenzo, Science, 270 255 (1995) 

[31] I.L. Chuang, R. Laflamme, PW. Shor, and W. H. Zurek, Science 270 1633 (1995) 
[32] W. G. Unruh, Phys, Rev. A 51 992 (1995). 

[33] R. Laflamme, C. Miguel, J. Pablo Paz, and W. H. Zurek, Phys Rev. Letters bf 77 198 
(1996). 

[34] R. Feynman, Int Jour. Theoret. Phys. 21 467 (1982). 



19 



[35] S. Lloyd, Science 273 1073 (1996); C. Zalka Efficient Simulation of Quantum Sysytems 

by Quantum Computers, Report No. BUTP-96/11, March 25, 1996, Los Alamos 

Archives |quant-ph/ 9603026] . 
[36] P. Shor, in Proceedings of the 35th Annual Symposium on the Foundations of Computer 

Science (IEEE Computer Society, Los Alamitos, CA 1994), p. 124. For a recent review 

see A. Ekert and R. Jozsa, Rev. Mod. Phys. 68 733 (1996). 
[37] P. W. Shor, Phys. Rev. A 52 R2493 (1990); R. LaFlamme, C. Miquel, J.P Paz, and 

W. H. Zurek, Phys. Rev. Letters 77 198 (1996). 
[38] S. Lloyd, Phys. Rev. Letters, bf 75 346 (1995); D. P. DiVencenzo, Phys. Rev. A 51 

1015 (1995); D. Deutsch, A. Barenco, and A. Ekert, Proc. Roy. Soc. London A 449 669 

(1995). 

[39] This sequence first appeared in the literature in conjunction with work on relaxation 
of hierarchical systems (B. A. Huberman and M. Kerzberg, J. Phys. A: Math. Gen. 18 
L331 (1985)) and recently in work on aperiodic Ising chains (F. Igloi and L. Turban, 
Los Alamos Archives |cond-mat / 9606118]) . The sequence is referred to as the hierarchical 
sequence. The author thanks Dr. Turban for these references. 

[40] R. Feynman, Optics News, 11 11 (1985); reprinted in Foundations of Physics, 16 507 
(1986). 

[41] P. R. Halmos, A Hilbert Space Problem Book, 2nd Edition, Springer Verlag, New York, 
1982. 

[42] M. Kolaf amd F. Nori, Phys. Rev. B 42 1062 (1990) 

[43] M. Queffelec, Substitution Dynamical Systems- Spectral Analysis, Lecture Notes in Math- 
ematics 1294, Editors A. Dold and B. Eckmann, Springer Verlag, Berlin 1987. 

[44] D. P. DiVicenzo and P. J. Steinhardt, Progress and Current Issues in Quasicrystals, in 
D. P. DiVincenzo and P. J. Steinhardt Eds., Quasicrystals the State of the Art Directions 
in Condensed Matter Physics- Vol. 11, World Scientific Singapore 1991, pp 1-15. 

[45] For example: M. Kohmoto, L. P. Kadanoff, and C. Tang, Phys. Rev. Letters 50 1870 
(1983); M. Kohmoto, International Jour. Modern Phys. 1 31 (1987); G. Gumbs and M. 
K. Ali, Phys. Rev. Letters 60 1081 (1988); M. Kolar and M. K. Ali, Phys. Rev B 39 
426 (1989). 



20 



FIGURE CAPTIONS 



Figure 1. Initial and Final States for Counting GQTM for the First 2 n Binary Numbers. 
All lattice qubits are in state |0) except those at sites and n + 1 which are in state |2). 
The initial and final head states are shown as wave packets with internal head states |0) and 
|1) to the left and right respectively. 



Figure 2. The Potential Distribution for N = Q. The potential height, in units of 
K{1 — 7), is plotted against path position. A potential height of 1 or 2 corresponds respec- 



tively to one or both terms of Eq. 11 being active. 



Figure 3. Energy Band Spectra for the Associated Periodic Systems for 7 = 0.999 for 
three values of N as a Function of Momentum k. The momentum range is from the least 
value for which \TrZ^\ ~ 2 up to k = arccos7. Within each long band or ribbon associated 
with each N, the upper horizontal line segments correspond to energy bands; the lower 
horizontal segments correspond to energy gaps. Very short bands or gaps appear as points. 
For N = 15 and especially iV = 18, band and gap density is too high in several momentum 
regions to be resolved at the scale shown. 

Figure 4. Energy Band Spectra for the Associated Periodic System for the Same Values 
of N as were used in Figure 3 but with 7 = 0.99. See the caption of Figure 3 for additional 
details. Note that, compared to the spectra in Figure 3, the gaps are much wider and the 
energy bands are narrower. 



Figure 5. Energy Band Spectra for the Associated Periodic System for 7 = 0.9 and 
N = 6, 8, 10. The caption of Figure 3 gives more details. Lower values of N are shown 
because at the higher values of N the fluctuations of TrZ^ are so extreme that calculations 
become impractical. Even for the values of N = 8 and iV = 10 there are no assurances 
that the grid used (with more than 7,000 equally spaced points) is fine enough to capture 
all bands. 



Figure 6. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N = 10 and 7 = 0.999. The momentum range is the same as that used in 
Figure 3. The corresponding energy band spectrum for the associated periodic system, from 
Figure 4, is shown at the bottom of the figure. 

Figure 7. An Expanded Plot of the Log (base 10) of the Landauer Resistance as a Function 
of k for N = 10 and 7 = 0.999 over the Range 0.034 < k < 0.040. The corresponding 
expanded band spectrum is also shown. The correspondence between minima in log LR and 
band-gap boundaries is clearly shown as is the decrease in the magnitude of the minima 
with increasing gap width. 

Figure 8. The Log (base 10) of the Landauer Resistance as a Function of the Momen- 
tum for N = 18 and 7 = 0.999. The corresponding energy band spectrum is also shown. In 
order to show values of log LR in the neighborhood of 0, the ordinate values were cut off at 



21 



20. Fluctuations in log LR are so rapid and extreme that individual maxima and minima 
in the neighborhood of are not resolved. 

Figure 9. Expanded Plot of the Log (base 10) of the Landauer Resistance for N = 18, 
7 = 0.999 over the range 0.032 < k < 0.033 The corresponding expanded band spectrum is 
also shown. The range is a small region in Figure 8 showing especially rapid fluctuations. 
The correlation between the height of maxima in log LR and widths of associated gaps is 
shown clearly. 

Figure 10. A Very High Resolution Plot of the Log (base 10) of the Landauer Resistance 
for TV = 18, 7 = 0.999 Over a Very Small Momentum Region of 0.03225 < k < 0.03228. 
The expanded band spectrum is also shown. Individual maxima and minima are clearly 
resolved. As is the case in Figure 7 for N = 10, the location of minima at or near band-gap 
boundaries, and the decrease in the minima with increasing gap widths are clearly shown. 

Figure 11. The Log (base 10) of the Landauer Resistance Plotted as a Function of the 
Momentum for N — 10 and 7 = 0.99. The corresponding band spectrum is also shown. The 
momentum range is from the least value of k for which \TrZ^\ ~ 2 up to arccos7. High 
maxima in log LR are clearly assocated with wide gaps. 

Figure 12. Plot of the Log (base 10) of the Landauer Resistance and Corresponding Band 
Spectrum for N = 8 and 7 = 0.9. The values of k range from the smallest value for which 
\TtZn\ ~ 2 to arccos7. The values of log Li? are cut off at 20 in order to show the fine 
structure at smaller values. 

Figure 13. Expanded Plot of the Log (base 10) of the Landauer Resistance and Corre- 
sponding Band Spectrum for N = 8 and 7 = 0.9 for 0.400 < k < 0.415. The association 
between minima in log LR and band-gap boundaries is weaker than for values of 7 closer to 
1 in that some minima are close to the middle of gaps. 
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